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Abstract. We present an exhaustive analysis of the light-induced frequency shifts 

ps^ I of the photoassociation lines of ultracold metastable '^He* atoms in a magnetic trap. 

K*" ' The measurements of the shifts of several vibrational levels bound in the purely long- 

^^ . range J = 1, 0+ potential linked to the 2^5'i — 2^Po asymptote were reported in 

»>f-\ ' a previous paper [J, and are analyzed here. The simplicity of this system makes 

^^ ■ it very appropriate for a detailed study. Indeed, the purely long-range character of 

>^ , the excited potential allows one to calculate exact excited molecular wavefunctions 

^— N ' and to use asymptotic expansions at large internuclear distances of the ground state 

*>>; . wavefunctions appearing in Franck-Condon type integrals. Moreover, the number of 

C^ I collisional channels to be considered is strongly reduced by the absence of hyperfine 

^ ' structure for '^He* and the use of polarized ultracold atoms and polarized light. This 

allows us to derive semi-analytical expressions for the shifts showing explicitly their 

r*^ ' linear dependences on the s-wave scattering length a of spin polarized metastable ^He* 

O ■ atoms. This explains how it is possible to derive the measurement of a from these shifts. 
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1. Introduction 

Since tlie breakthrougii of tlie methods for cooling and trapping atoms twenty years ago, 
a renewal of interest occurred for photoassociation (PA) experiments, in which a pair of 
colliding cold atoms is promoted into an excited molecular bound state by absorption 
of laser light. Research in this domain provided a large amount of high precision results 
in the field of molecular spectroscopy and cold atom collisions J2I • When the intensity 
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of the PA laser light is high enough, it is likely to induce frequency shifts which have to 
be controlled for accurate spectroscopy. Such frequency shifts are interesting processes 
which result from the perturbation of molecular states by the electromagnetic radiation. 
They have been studied by several groups dealing with cold alkali atom clouds, both 
experimentally and theoretically. In particular the cases of Na |31 E] and Li E] 
have been investigated, and the theoretical background has been provided by [71 |H]. 
Even if the physics behind these experiments is well understood, a detailed comparison 
between theory and experiment is somehow limited by the incomplete knowledge and 
the complexity of the molecular potentials for alkalis. 

In the present article we investigate theoretically the light-induced shifts in the 
photoassociation of ultracold 2^5*1 metastable ^He* atoms. This case is significantly 
simpler than the case of alkalis for several reasons. First, '^He* has no hyperfine structure 
so that the number of channels involved in the calculation is smaller. Second, we 
consider here photoassociation processes which produce excited molecules in the purely 
long-range Oj potential which can be exactly expressed in terms of atomic parameters 
(atomic C3 dipole-dipole coefficient and fine structure of the 2^P excited atomic state). 
The wavefunctions of the excited molecular states which describe giant dimers can 
thus be exactly calculated CHI El- Third, since the wavefunction of the excited 
molecular state is non-zero only for large values of the internuclear distance r, one can 
use asymptotic expansions of the ground state molecular wavefunctions in the Franck- 
Condon type matrix elements which appear in the expression of the shifts. We thus do 
not need to know the wavefunction in the short range domain where it depends critically 
on the ground state potential. We show in this article that it is then possible to obtain 
semi-analytical expressions for the light-induced shifts showing explicitly their linear 
dependence on the s-wave scattering length a which describes the collisions between 
spin polarized ultracold metastable ^if e* atoms. 

The results derived in this paper were actually used for the interpretation of a 
PA experiment with ultracold metastable '^He* atoms, which was recently published 
by our group P], and from which we derived a value for a. Actually, such light shift 
measurements were a first step to finding a reliable value for a, an important parameter 
ruling the interatomic interactions and the Bose Einstein condensate properties. The 
derivation of a from the light shifts was a preliminary study to reaching an even much 
more accurate value by a second method ^2j- This was based on two-photon PA 
spectroscopy and dark resonances, from which we measured the binding energy of the 
least bound state in the interaction potential of the two metastable atoms. It allowed 
us to find a value for a at least a hundred times more precise than the best previous 
determinations, in disagreement with some of them. The preliminary one-photon PA 
measurement of a, derived from the understanding of the light shifts, was very helpful to 
focus the frequency scan for the final two-photon experiment to a restricted range centred 
around a reliable value. The interpretation of the physical origin of the light-induced 
frequency shifts has been briefiy explained in the article reporting our experimental 
results [T]. The detailed analysis given here is exhaustive and justifies the derivation of 
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the value of a from our previous measurements. 

We first explain in section |21 the physical origin of the shift of the PA line. We 
show that it can be related to the well known problem of a bound state coupled to a 
continuum of states and to a set of bound states [12] . In our problem of ultracold atoms 
held in a magnetic trap, the unique significant entrance coUisional channel leading to real 
photoassociation processes corresponds to two spin polarized atoms colliding in an s- 
wave {£ = 0). We point out the two possible origins of the dependence of the shift of the 
PA line on the scattering length a characterizing the collision in an s-wave: dependence 
on a of the binding energy of the least bound state in the potential of the entrance 
s-wave channel and dependence on a of the Franck-Condon overlaps of the excited and 
ground molecular wave functions. We emphasize also the importance of taking into 
account other non s-wave ground state channels for evaluating the shift of the PA line. 
This is due to the fact that the shift is sensitive to non resonant couplings to these 
channels which may have a non negligible contribution, whereas the photoassociation 
rate is due to resonant couplings occurring only in the s-wave entrance channel. Since 
the i ^ ground channels coupled to the excited bound molecular state must satisfy 
polarization selection rules, one expects a laser polarization dependence of the shifts. 

We apply these considerations in section 01 to the photoassociation of two spin 
polarized ^He* atoms in a well defined rovibrational state of the excited J = 1,0^ purely 
long range potential. Using the selection rules of electric-dipole transitions, we identify 
the various ground state channels to be considered for two different laser polarizations. 
We also give the expansion of the excited molecular state in the Hund's case (a) basis 
which is appropriate to the calculation of the matrix elements of the molecule-laser 
coupling . 

The explicit expression of the shift is derived in section E] from an effective 
Hamiltonian approach. Two formulas previously obtained by different methods are 
shown here to be equivalent. In the first one, the shift is given by the well known 
second-order perturbation expression jTSl E] • The second one is easier to use since it 
involves only the wavefunction of the excited state and ground state wave functions at 
a given scattering energy with well defined boundary conditions [HI |T3l Uni Hj. The 
fact that the excited molecular states in the potential 0^ are long range allows us to 
use the asymptotic form of the ground state wavefunction. The second expression of 
the photoassociative shift derived in section |3] can thus be simplified and in section j^ we 
derive expressions for the shifts 6"" of three PA lines (f = 0, 1, 2) showing explicit linear 
dependence of 6"" on the scattering length a. The curves giving 6^^^/ 6""^^ and 6^^'^/6^^'^ 
as a function of a are given. These theoretical curves were used in ^ to deduce a range 
of values of a from our measurements of the shift ratios. 
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Figure 1. Bound state \(j)b) of energy eb coupled to a set of bound states {(pi) of energy 
Ei, i — 1 . . .n and to a continuum of states \ipE') of energy < E'' < +cx). Ttie matrix 
element which couples |0b) and a state |'0B') of the continuum is noted Qi,{E'), and the 
one which couples 10;,) and a bound state \4>i) is noted Clf,i. These couplings induce a 
shift of the energy Cb of the excited state (pb- The continuum of states of energy £" > Cb 
gives a contribution to the shift towards negative energies, whereas the continuum of 
states of energy E' < eb and the bound states give a contribution to the shift towards 
positive energies (see equations ^ and ^). 



2. Physical origin of the shift 



2.1. Shift of a bound state coupled to a continuum and to a set of bound states 

We first recapitulate here on a few results about the coupling of a bound state to a 
continuum and to a set of bound states [121 ■ Let 0b be a bound state of energy e^,, 0j 
a set of n bound states of energy Ei < 0, i = 1 . . .n, and ipE' a continuum of energy 
normalized states of energy < E' < +oo (see FiglH). We assume that the states are 
chosen so that the Hamiltonian H is diagonal in the two subspaces spanned by (pb on 
the one hand, and the (/>i and ipE' on the other hand : 

{(Pb\H\^E') =^b{E') 
{(t>b\H\(j)i) = Qbi 
{(j)i\H\(j)j) = Ei6ij 
{iIje'\H\^e") = E'6{E' - E") 
{(j)i\H\^E') =0 
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We have taken h = 1. We note for further convenience T{E') = \Qf){E')\'^. The rate 
r(£") has the dimension of an energy as the continuum states ^pE' are energy normahzed. 
The probabihty of excitation of the bound state 06 from a continuum state ipE exhibits a 
resonance when E = eb + Si,. The shift Sf, = ^^°"-* + ^tound ^g ^.j^g g^^^ ^^ ^.j^g contributions 
^cont q£ ^Yie continuum states "^e' and 5^°""°' of the bound states (pi. 

^cont ^-py f ^E' ^^'^l (1) 

Jo ^h - E' 

^.ouna ^ J2 ^ (2) 

t 

In (P) PV denotes Cauchy's principal value integral. We find that that the continuum 
of states of energy E' > eb gives a contribution to the shift towards negative energies, 
whereas the continuum of states of energy E' < tb and the bound states cpi give a 
contribution to the shift towards positive energies. 

2.2. Application to photoassociation 

We now apply the previous considerations to the photoassociation (PA) of a pair of 
colliding atoms in a molecular bound state of energy Eb in an excited potential. In 
the light-dressed picture, the excited potential is shifted by an amount —ui (we remind 
that h = 1) corresponding to the energy of a photon of the PA laser. The problem is 
analogous to the one of a bound state coupled to a continuum and to a set of bound states 
(see Fig. El in the metastable helium case). Using the notations of paragraph 12. H cpb is 
the excited bound state of energy eb = Eb — uj, the states ipE' and ipi correspond to the 
scattering and bound states in the ground state potential, Qb{E') and Qbi are the light- 
induced couplings between the ground and excited states, which are proportional to the 
square-root of the PA laser intensity. We assume that these couplings are weak enough 
so that (pb can be considered isolated from the other states in the excited potential. We 
also assume that only the least bound state of energy En in the ground state potential, 
which is the closest to resonance, has a non negligible influence on the spectroscopy 
of the excited state. In the context of ultracold collisions between atoms, the energy 
of the initially populated scattering state E^o is close to the dissociation threshold of 
their interaction potential, which means E^o — 0. The unperturbed resonance condition 
Eoo = £(, imply that tb — 0. Using these simplifications, the contributions in equations 
dH) and (J2I) give for the total shift: 

Sb^-Vvj^^dE'^-^-"^^ (3) 

As T{E') and Vt^^ are proportional to the laser intensity, the shift 5b is also 
proportional to the laser intensity. 
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Figure 2. (a) States and potential curves involved in the experiment of Pj. Since 
we are in the dressed-atom picture, the difference between the dissociation energies of 
the two potential curves is equal to the detuning A = oj — ujq, where lo JI-k is the laser 
frequency and loqIItx is the frequency of the atomic transition. The spin polarized 
""He* atoms collide with an energy E^o in the ^E+ potential, which supports 15 bound 
states. As the scattering length a for this potential is large and positive, the energy 
of the least bound state -y = 14 is close to the dissociation threshold. The atoms are 
photoassociated in bound states u = 0, 1, 2 of the purely long-range J ~ 1,0+ potential. 
In this figure, the laser frequency lo jlix is close to resonance with the transition which 
drives the colliding atoms to the w = state, (b) Problem in terms of bound states 
coupled to a continuum and a set of bound states. In the case where the coupling 
between the ground and excited states is small compared to the level spacing in the 
excited potential, the influence oi v — \ and w = 2 on the photoassociation of u = 
can be neglected. This scheme is then analogous to that described in Figure ^ Here 
-y = is the bound state |0f,) which is optically coupled to the scattering states \'4i'^ 
and to the bound state \4>i^\^ of the ^I]+ potential. 



2.3. Coupling with a ground state s-wave channel 

Considering the coupling with an s-wave coUisional channel characterized by a scattering 
length a, two origins can be identified in equation Q for the dependence of the shift on 
a. 

First, in case of a large and positive scattering length a, the binding energy En of 
the least bound state is approximated by — 1/2/ia^, where /x is the reduced mass of the 
colliding system. As noticed in |8^, and also mentioned above, the contributions from 
the bound state and from the continuum of states have opposite signs in equation 0: at 
zero scattering energy, the continuum contribution is negative, whereas the bound state 
one is positive, and all the more important that the scattering length a is large. The 
shift is sensitive to the vicinity of the least-bound state in the ground state interaction 
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Figure 3. Wavefunctions w = 14 for the molecule in the ground state potential ^S+ 
and w = 0, 1, 2 for the molecule in the excited purely long-range J = 1,0+ potential. 
The overlap of the wavefunctions w = 14 with v = 0, 1, 2 appears in ^1,(1=14:) in the 
bound states contribution (O to the shifts of 1; = 0,1,2. The wavefunction of 1; = 1 
crosses zero an odd number of times whereas the wavefunctions for v = and v = 2 
cross zero an even number of times. The wavefunction oi v = lA varies monotonously 
in the range of interatomic distances where the wavefunctions of v = 0, 1, 2 take non 
negligible values. As a consequence, in case of the overlap with u = 1, positive and 
negative contributions approximately cancel each other, and the contribution of w = 14 
to the shift of w = 1 is expected to be smaller than its contribution to the shift oiv — 
and V = 2. 



potential and thus to a. 

The second origin for the dependence of the shift on a is found in the stimulated 
rates f2^^ and T{E') which contain overlap integrals involving a-dependent wavefunctions 
for the electronic ground-state and wavefunctions for the excited molecular state. The 
shift for each excited molecular state depends differently on a. Figure El shows for 
example the wavefunctions of the ground and excited molecular bound states involved 
in the experiment described in |1 . Due to the fact that the molecules in the excited 
electronic state lie in a purely long range potential, the nodes of their wavefunctions are 
localized at a distance corresponding to the vanishing exponential tail of the molecular 
wavefunction in the electronic ground state. As mentioned in jB|, if the wavefunction of 
the excited bound state changes sign an odd number of times, positive and negative 
contributions approximately cancel each other in the bound-bound overlap integral 
contained in fif,^ , and the contribution Q to the shift is smaller than in the case of an 
even number of changes of sign of the wavefunction. Therefore, the investigation of the 
shift for several excited states provides complementary information on the scattering 
length. 
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2.4- Contribution of other ground state channels 

At very low scattering energy E^o, atoms colliding in a £ 7^ partial wave are held apart 
by the centrifugal barrier i{i + l)/2/ir^, which prevents them from being efficiently 
photoassociated. The free-bound stimulated rate TilEoo) > for an incoming i > 
partial wave can be neglected compared to rf=o(-Eoo)- Only s-wave coUisional channels 
contribute to resonant photassociation processes producing excited molecules. However, 
the excited level shift (jHI) involves the contribution of the whole spectrum. States 
of higher energy E', for which Ti^qIE') is not negligible compared to Ti^qIE'), can 
contribute noticeably to the shift even though they are non resonant. For example, in a 
previous experiment with a Bose-Einstein condensate of Na [H] , the d-wave contribution 
is enhanced by a shape resonance and is much larger than the s-wave contribution. 

Owing to the contribution of non s-wave collisions, the shift is expected to depend 
on the polarization of light, as already mentioned in ^7]. The collisional channels which 
are coupled to the excited state depend on the polarization of light. Therefore the 
investigation of the shift as a function of the light polarization allows one to tune the 
relative contribution of the s-wave collisional channel, which depends on a, and of the 
other i ^ contributions. 

As a consequence, when investigating the light shift of the excited state, one has to 
carefully take into account all the ground state collisional channels which are coupled, 
and not only those which are resonant. 

3. Helium interaction potentials 

3.1. Collisions between metastable helium atoms 

In the absence of a light field and neglecting the weak spin-spin interaction between 
atoms, the total electronic spin 5* and the relative angular momentum i of two spin-1 
colliding '^He* atoms are conserved to good approximation. A basis to describe the 
collision of a pair of colliding atoms a and b is given by (JD), where Sa, Sb are the 
spins of atom a, atom b, i the relative rotational angular momentum, Ms and M^ are 
the projections of S and i on the lab-fixed quantization axis. They have no electronic 
orbital angular momentum as the atoms are in a 2^S state. 

\{Sa = hSt = l)S,Ms;e,Me) (4) 

As shown in [Appendix A[ due to the bosonic nature of the metastable helium atoms, 
singlet {S = 0) and quintet {S = 2) spin states colhde in even £ partial waves, whereas 
triplet {S = 1) spin states collide in odd i partial waves. Due to the bosonic nature 
of both the helium atom and helium nucleus, singlet and quintet (triplet) states have 
the gerade (resp. ungerade) symmetry under the inversion of the electronic coordinates 
through the center of charge of the molecule. As both atomic 2^5"! orbitals are spherical, 
the spatial part of the electronic wavefunction is not changed by a symmetry through 
a plane containing the intermolecular axis, and the colliding pair has a + symmetry. 
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Figure 4. Potential energy curves involved in the collision between metastable helium 
atoms as calculated bv [THl ITOl 0011^1 1^ . 



The colliding states are E states because the colliding metastable helium atoms have 
no electronic orbital angular momentum. Figure |3] shows the interaction potentials 
^S+,'^E+ and ^E+ associated to the quintet, triplet and singlet states respectively. 

3.2. Ground states optically coupled 

In the experiment described in PJ, the ^ife* atoms are spin-polarized and held in a 
magnetic trap. The atoms collide at energies such that only the s-wave part of the 
scattering wave contributes to the photoassociation, and the state which characterizes 
the internal and rotational parts of the system is \{Sa = I, Sb = I), S = 2, Ms = 2;i = 
0,M£ = 0). The total angular momentum J is defined as S + L + i, where L is the 
total electronic orbital angular momentum and L = for a pair of metastable atoms. 
Therefore the colliding state is a J = 2, Mj = 2 state. 

The total angular momentum J is not conserved in presence of a laser field. We 
note J' and Mj/, the total angular momentum of the excited state and its projection 
on the lab-fixed quantization axis. We remind here the selection rules for electric-dipole 
transitions : 

• |^-i| <J'< \J + M 

• Mjr — Mj = q where q = —1, 0, 1 for cr~,7r and a~^ polarization respectively 

• there is no coupling between states of the same symmetry under the inversion of 
the electronic coordinates through the center of charge of the quasi molecule 

In the experiment described in pQ, the atoms are excited in vibrational levels of 
the purely long-range J' = l,Mj/ = 1;0^ potential. As already mentioned at the end 
of section El the initial scattering state \{Sa = I, Sb = 1), S* = 2, Ms = 2;i = 0, Mi = 0) 
is not the only ground state which will be relevant for evaluating the light shift. Once 
the excited 0^ state is formed, it can be optically coupled to other ground states which 
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contribute to the light shift. Because of the selection rules, the projection of the total 
angular momentum in the ground state is Mj = Ms + Mi and has to be Mj = 1 — q. 
Furthermore, S and i should be such that they can form a total angular momentum J 
that satisfies the triangular inequalities with J' = 1 and the photon angular momentum 
equal to 1. Therefore, the properly symmetrized states (see (jA.lj) in [Appendix A| ) which 
are coupled to the excited state by a^ and a^ light are given by (0) and (jH)) respectively 
using the simplified notation \S,Ms]i,Mi) (see f)A.2|) ): 

(5) 



|00;20) (6) 



Several remarks can be formulated : 

• In ® M5 + M^ = 2 and in (0) M5 + M^ = 0. 

• The scattering wavefunctions for the channels |22;00) and |20;00) which appear 
in both © and (jH)) depend on the quintet scattering length a which is searched 
in pp. The wavefunctions in the other channels do not depend on a, because they 
correspond to either a singlet state, or to d- or (7-wave collisions. 

• By changing the polarization of light, one can change the relative contribution of 
a-dependent and a-independent channels as the coupled channels are different in 
both cases. 

• The state |00; 00) involves the scattering length of the ^S+ potential (see Fig. HJ, 
which is estimated to be ai = 35± 15 a^ according to previous calculations [TH1 I23J. 

3.3. Description of the excited state 

Contrary to the case of molecular potentials connected to the S + S asymptote 
investigated in paragraph 13. 2^ S and i are no longer good quantum numbers to describe 
the molecular potentials which dissociate into the S + P asymptote. Due to the fine 
structure splitting of the P state, the spin S is mixed to the orbital angular momentum 
L and is not conserved. The interaction between the S and P atoms is not isotropic as 
it depends on the relative orientation of the electronic orbital angular momentum and 
the molecular axis, and i is not conserved. The remaining good quantum numbers to 
characterize the excited state are J,Mj and Q, the projection of L + S* on the molecular 
axis. The purely long-range 0^ potential used in [l] was previously investigated in great 
detail experimentally and theoretically [01 dUl E] . We recapitulate here the essential 
points needed for our purpose. 

The purely long-range 0^ potential which connects to the 2^5*1 + 2^Po asymptote 
is mainly determined by the atomic C3 dipole-dipole coefficient and the fine-structure 
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splitting of the 2^P state. It is shallow (about 2 GHz deep) and supports five bound 
states which have inner turning points at internuclear separations of approximately 
150ao. Their electronic state is close to that of two separated 2^5" and 2^Po atoms. The 
decomposition of the electronic 0^ state in the Hund's case (a) basis is needed to estimate 
the light-induced coupling between the ground and excited states. It has contributions 
from the states ^Ilu,^Ilu,^'^t ^^^ ^S+. The ^S+ state is not present as it belongs to the 
0~ subspace. We write the decomposition of the purely long-range |0^) electronic state 
in the Hund's case (a) basis in equation ((Zj), where the r-dependent a coefficients are 
obtained by the diagonalisation of the r-dependent electronic Hamiltonian ^T] : 

K) = «Bs+(r)|^S+) + asu^{r)\'U^) + mn^WI^H^) + a.^+WI^S^) (7) 

Experimentally, the rovibrational J = 1 states are probed, where J is the total 
angular momentum of the molecule. The J = 2 states do not exist as they are not 
symmetric under the exchange of the two bosonic helium nuclei, and the J = 3 states 
could not be observed experimentally because of the very weak coupling with the initial 
state. This can be understood by the fact that the transition is forbidden in the limit 
of two non-interacting atoms and that the electronic state of the excited long range 
molecule is close to the one of two non interacting atoms. As the J = 1 states are 
experimentally produced from a J = 2 state (see subsection I3.2|l . owing to selection 
rules, the projection Mj of J on the molecular axis is Mj = 1. Due to rotational 
couplings, more Hund's case (a) states than the four previous one are involved in 
the decomposition of the photoassociated levels, but their contribution is only non 
negligible for interatomic distances where the purely long-range 0^ potential crosses 
other adiabatic potentials. This occurs at interatomic distances well below the inner 
turning point of the photoassociated levels. Therefore, we keep only the four electronic 
states ((Zj) in the expression of the wavefunction for the rotating molecule. 

\J' = l^Mj. = 1;0+) = S,a,(r)|J' = l,Mj. = l,i) (8) 

where i = ^S^,^ H„,^ H„,^ SjJ^. In Appendix B[ we give the explicit expression ()B.4|) 



of the states \J' = l,Mj' = l,i) in terms of states of the Hund's case (a) basis for a 
rotating molecule | J, Mj, fl, S,T,, A, w) . 

4. Two equivalent formulas for the molecular light-shift 

We derive here formulas for the shifts of the bound levels in the electronic excited state 
I J' = l,Mj' = l;0jj^) described in section EISl due to the optical coupling with ground 
state channels described in section 13.21 Using an effective Hamiltonian treatment and 
two different expressions for the field-free propagator for the ground state channels, we 
show the equivalence between two formulas established in previous papers. The first 
one involves a sum over the whole energy spectrum in the electronic ground state and 
was derived in ^31 and in [Hj in the context of photoassociation. The second one derived 
by different means in ^J^J^lIT] involves only ground state wavefunctions at a given 
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scattering energy, and is therefore much easier to estimate numerically. The equivalence 
between the two formulas which we demonstrate here can be regarded as an extension 
of the sum rule described in J21I. Let us emphasize that the contribution of the bound 
states to the shift, which is explicit in the first formula, is also included in the second 
one. 

We now present the effective Hamiltonian treatment used to derive the shifts. 
Given that the couplings between the ground state channels due to the very weak spin- 
dipole interaction |2S1 can be neglected, we restrict ourselves to a two-channel treatment 
between a \g) = \S,Ms;i,M() ground state and the |e) = \J' = l,Mj/ = 1;0^) excited 
state. This approximation will be discussed at the end of this section. The shifts 
Ss,Ms;i,Mt calculated in this way for the states (0),® can then be added up to calculate 
the total shift S: 

S= J2 ^^ (9) 

g=S,Ms-AMe 

The coupling induced by the laser of intensity / between a ground channel g = 
{S,Ms;i, Mi) and the excited channel e is defined as ^s,Ms;i,Me{r) = f2eg(r) and 
is calculated in appendix B. Note that ^eg{f) is obtained after an integration over 
electronic variables and over the angular variables of the molecular axis so that fieg(^) 
acts only on the variable r, distance between the two nuclei. The problem can be 
described by a r dependent two-channel Hamiltonian H2C ^^ '■ 

H2C = H^ + Vr (10) 

H - 1 * " + I ^^^") ° 



V, 



R 



2/irfr2 


+ 
\ 





fige( 


^eg{r) 






In flTUj) Vg{r) is the ground state potential which is the sum of the interaction potential 
2S+i^+ ^^^ Qf ^j^g centrifugal potential £(£ + l)/2iir'^. Ve{r) is the J = 1, 0+ potential. 
1 is the identity matrix. 7 accounts for spontaneous emission from the excited state. 
We choose the reference of energy at the dissociation energy of the potential associated 
to the channel g : {S = 2, Ms = 2;£ = 0, M^ = 0), which corresponds to the collision 
between spin polarized atoms. A is the detuning u — ujq between the laser and the 
atomic frequencies (see Figure |21^a)). In presence of a magnetic field B, the potentials 
corresponding to channels with different projections Ms are shifted proportionally to 
gsfJ'sB, where gs is the electron spin gyromagnetic ratio and ^b the Bohr magneton. 
We define therefore the energy E"^ and the wavevector fc^ which will appear in the 
matrix elements of the effective Hamiltonian Heff{E) : 

Ems = ^s(^) = (Ms - 2)gsfiBB (11) 



kZ, = \ME-E^) (12) 
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To obtain the effective Hainiltonian, we define P and Q as the projectors onto the 
states in the ground channel and the excited channel respectively p3|) . These states can 
be characterized by the quantum numbers a and (3 which denote either a vibrational 
number or the scattering energy of the eigenstates of PHqP and QHqQ. 

p = X^ifl',«)(fl',«i 

a 

Q= ^|e,/3)(e,/3| (13) 

Defining the field free propagator at energy E with outgoing boundary conditions 
Gq{E) = {E — Hq -\- ie)^^ (chapter XIX in ^^), we obtain the effective Hamiltonian 
Hi,ff for the excited state: 

HeffiE) = -^^ + K(r) - (A + ^|) + QVnPGi;\E)PVnQ (14) 

From now on, we will consider that only one bound state of vibrational number 
V of the excited potential is efficiently coupled by light, which corresponds to the 
isolated level approximation. Its validity will be discussed later. We have therefore 
Q = Qv = |e, f)(e, f |. The matrix element {e,v\Heff\e,v) gives the energy of the bound 
state V perturbed by light. The following expression which appears in {e,v\Heff\e,v) 
contains a real part S^{E) which corresponds to the light-induced energy shift of the 
level V, and an imaginary part r'"{E) corresponding to the light-induced broadening: 

{e,v\Q,VRPGi-'\E)PVRQ,\e,v) = 6^g{E) - iTl{E)/2 (15) 

Let us consider the matrix element {t'IPGq (E)P\r) = Q^(E]r,r') of the propagator 
PGq {E)P. It is related to the matrix element of (fT5|) by equation (|T6|l . noting 0^(r) 
for the wavefunction associated to the state |e, f) which is real. We use also the fact 
that Vteg is real so that Vteg = ^ge- 

/•OO /"OO 

{ci>l\QVRPGf\E)PVnQ\ct>l) = / n,g{r)n,,{r')g+{E-r,r')<Pl{r)ct>l{r')drdr' (16) 



g^{E; r, r') can be written in two following ways : 
(i) The definition Gq (E) = [E + it — Hq)^^ of the propagator and the expression 



^(+) 



of PGq {E)P in terms of discrete eigenstates \g,i = 1 . . .n) of energy E^, and 



scattering eigenstates \g,E') of energy E' in the channel g provide (fTTjl . In (fTTjl . 
VV denotes Cauchy's principal value, 0g(r) is the real wavefunction associated to 
the bound level \g,i = 1 . . .n), 0^(r) is the energy normalized real wavefunction 
associated to the scattering state \g, a = E) and E^ the dissociation energy of the 
potential associated with the channel g (see pi|) ). 

g+{E;r,r) = 2_^ E - E^ / oo E - E' — ^^ ~ ^^'^s ^'^ ^'^s ^ (17) 
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The wavefunction (pf is the solution of the radial Schrodinger equation associated 



with the channel g, and has the following boundary condition for r 
involves the phase shift Tig{k'^) and the wavevector k'^ 



oo, which 
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OO 
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■nkf 



sm 



TX 



k^r-i- + r^gik 
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(ii) The solution of the differential equation which defines the propagator Q^{E]r,r') 
with proper boundary conditions gives (chapter XIX in [2S|) : 

^+(E;r,r') = -vr0f(r<)0f+(r>) (19) 

r< and r> are the smaller and the larger of r and r'. 0f'"'" is the solution of the same 



Schrodinger equation as (j)^, but with outgoing boundary conditions at infinity: 



kE+i 



CX3 



2/i 



exp 



vr 



kTr-ii^ + Vaik 



(20) 



Inserting ()17p and (fT^ into ()16|). and identifying the real and imaginary part of 
(fTHj) with the ones of (jT^. we get the light-induced broadening and shift. Both (jTTj) 
and (fT^ give the same form for the light-induced broadening ^"{E): 

2 



riiE) = 27r 



0f(r)fie,(r)C(r)rfr 



(21) 



The expression for the shift 6^{E) is different whether equation ()17p or p9|) is used, 
giving respectively (j^ and (J2S1): 



^m = y: 



r <Pl{r)n,,{r)<j>:{r)dr\ 



j=i 



E-E: 



vv 



K{E) = 2n 



9 

oo 



r0f(r)fie,(r)0:(r)rfr 



(ir 



c?rU,(r)fie,(r')0f(r)0f(r')0:(r)0:(r') 



-t/E' 



(22) 



(23) 



Here, (pf is the energy normalized solution of the radial Schrodinger equation at 



energy E for the channel g, which is related to the real part of (j2(jp : 

0|;(r ^ oo) = 



^^ cos (k'^r-£^ + ri,ik: 



nk'^ 



The expression ()23|) , which was derived by different ways in f!M IT^ ITHl [7] is easier 
to use than the equation (j22j) as it involves only wavefunctions at a given scattering 
energy E rather than at all energies E' and E^ . Both formulas are derived from two 
exact expressions of the field-free propagator (fT7|) and (jT^ , and are thus equivalent j21] . 
In particular, (j^Hj) includes the contributions of the bound states to the shift, which are 
explicit in (1^ . 

Two approximations were done here. First, the isolated level approximation : 
when the PA laser is close to resonance with an excited bound level, the influence 
of the other excited levels is neglected. It is valid as far as the laser intensity is such 
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that 6'"{E) and r^{E) are small compared to the excited level spacing ^J^jn^l- If 
this approximation breaks down, the level shift is no more a quadratic function of the 
couplings fleg{r), and thus no more proportional to the laser intensity jHj. We are 
only interested here in shifts linear with laser intensity as observed in the experiment 
described in fD. Second, we neglected the spin-dipole interaction between ground state 
channels. Taking them into account would involve operators Q^VrP'Gq {E)PVjiQy 
in the effective Hamiltonian ^^, where P' is defined similarly as in p3|) but with a 
different ground state channel g'. The matrix elements of these operators would have to 
be included in the shift calculation. Their importance compared to the non neglected 
matrix elements of QvVrPGq {E)PVrQv is related to the probability of spin relaxation 
compared to the probability of elastic collision per collision event. These probabilities 
were calculated in jSHl 123 EH] and indicate that the matrix elements we neglected are 
four orders of magnitude smaller than those involved in our calculation of the shift. 

As a final comment, we notice that the Hamiltonian ()1U|) describes an optical 
Feschbach resonance JH]- The two potentials are shifted from each other using a laser 
field, similarly to what is achieved with a magnetic field to tune the scattering lengths 
pn] . The use of light coupling to tune the scattering length has already been investigated 
both theoretically [IHlliniEni and experimentally (SUES- The case of metastable helium 
was considered in [33]. Here, the shift is expressed as a function of the unperturbed 
wavef unctions, and will be evaluated as a function of the unperturbed scattering length 
in section |31 

5. Calculations of the shift for Hehum 

The expression of the light-induced energy shift in equations (jHI) and ()23|1 involves the 
overlap of the regular and irregular scattering wavefunctions with the wavefunction of the 
excited bound levels. As these excited bound levels lie in a purely long-range potential, 
the integrals do not depend on the short-range part of the ground state wavefunctions, 
and numerical calculations show that the integrals can be cut at Rmin = 100 Oq and 
Rmax = 500, 800, 1000 qq for v=0, 1 and 2 (see Figures El and Et)- Asymptoptic 
expansions for the ground state wave functions can therefore be used. Their accuracies 
were checked by comparing the results obtained in this way with the shift calculated 
using the wavefunctions of ab initio potentials [IHl EUl 1211 I22j- The calculation of the 
couplings ^sMs'AMti''^) ^^ reported in appendix B. 

5.1. Contributions of s -wave ground states 

In the experiment described in |Tj the excited state is probed by the excitation of atoms 
colliding in the S = 2, Ms = 2;i = 0, Mi = channel with a scattering energy E^o using 
cr~ light. The contribution of this channel to the shift involves the regular and irregular 
solutions of the radial Schrodinger equation at energy E = E^o = (^m5=2)^/2/^ for the 
^Si" potential and ^ = 0. These solutions are shown in Figure Eb- In the experiment 
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Figure 5. Wavefunctions involved in the calculation of the shift (|2Sl- Inset (a) 
shows the wavefunction of the excited level u = in the J = 1,0+ potential which 
takes non negligible values between Rmin = 100 oq and Rmax = 600 oq. Inset (b) 
shows the regular and irregular wavefunctions iJj2.2-q,o and ^"2,2:0.0 for the collision 
in an s-wave scattering channel at low scattering energy Eoo — (fcMs=2)^/2/^- The 

mm <r< Rmax- As a consequeucc, 

a) and 



^Ms=2\ 



wavevector k'^^^^ i^ ^^^h that /cms=2'' << Hot: R 

for this range of values of r ^^2, 2:0,0(7") « sin (fc^^^2('' ~ o)) 

''/'2.2;0,o('') oc — cos (fc^ 2(^ ~ '^)) — ^1- I^sct (c) shows the regular and irregular 

wavefunctions '4'S,Ms;e=2;M( and ^s,Ss;i=2,Mt fo^ ^^^ collision in a d-wave scattering 

channel. In the range Rmin ^ f < Rmaxj the interaction between the atoms colliding 

in a ^ wave is mainly due to the centrifugal potential £{£ + l)/2/ir^. 'ips,Ms-/=2,Me 

and ips,Ms;t=2,Mi are well approximated using spherical Bessel functions of the first 

and second kind (see equations H28|l and l|29|l ). It can be noticed that the product of 

the regular and irregular wavefunctions which appears in H2^^(l is of the same order of 

magnitude in (b) and (c). This shows that the d-wave contribution to the shift may 

not be negligible compared to the s-wave one. 
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described in |lj, the colliding energy is small enough for having the s-wave phase-shift 

r75=2,Ms=2;£=o,M,=o determined by a {r]s=2,Ms=2;i=o,Me=o - -^Ms=2«) and A;^g=2^ << 1 
for Rmin < ^ < Rmax- We Can thus approximate the regular and irregular wavefunctions 
for the atoms interacting in the ^E+ potential by their asymptotic expansion : 
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^2,2;0,o('^) 



2/i 



'^^Ms=2 



sm 



[k 



Ms =2 



r — a 



2/x 






COS [kt,„_o[r - a 



.00 / 
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Using equation (jSHI) and developing (j^ and (jS^I) for k 
equations : 



00 

Ms=2 



(24) 

(25) 
gives the following 



"2,2:0 



,o(-^o 



-4/i / dn / f/r2fi2,2;o,o(ri)0r(ri)fi2,2;o,o(r2)0r(r2) 



1^2 - a) 



1 + 



(^Ms=2)%,^ „,2 , (^2 -a)' 



■((ri-a)^ + 



(26) 



Equation ()26|) shows that when /c^ =2 ~^ O5 the s-wave contribution to the shift contains 
terms linear in a and quadratic in fc^ =2 ^"^^ thus linear in E^o- In the limit k^^ 2 = 0, 
equation (jTTj) shows explicitly the linear dependence of ^g 2*0 as a function of a: 



-4/i / dri / dr2n2,2-flfi{ri)(t)T\n)^2,2-flfi{r2)<PT\r2)r2 



^2,1:0,0(0) 



-4/i / dn / rfr2fi2,2;o,o(ri)0r(ri)fi2,2;o,o(r2)0r(^2) (27) 

Jo Jo 

The state jS" = 2, M5 = 0; £ = 0, M^ = 0) is also coupled by a^ light to the excited 

state. Its contribution to the shift involves wavefunctions which have the same form 

as (j^ and pHjl . where ^^3=2 has to be replaced by ^^3=0 = ■\/2/i(i?oo + '^Qs^^bB) 

according to (fTT|) and (fT^ . In case of a large magnetic field, the approximation 



^5=2,Ms=2;^=0,Mf=0 



-k^g^^a and the inequality k^^^^^r << 1 for R„ 



< r < Rr, 



may no longer be valid, and numerical calculations of the wavefunctions should be used 
to replace the formula ()26|1 . The same remark holds for the state IS* = 0, Ms = 0]i = 
0,Mi = 0). Additionally, the phase shift in this channel is related to the scattering 
length ai which is estimated to be 20 aQ < ai < 50 a^ according to previous calculations 
fHl I2n] • The uncertainty in ai induces an uncertainty on the calculation of the shift if 
a~^ polarization is present in the PA laser. 



5.2. Contributions of i ^ waves ground states 

The regular and irregular wavefunctions for a channel \S,Ms;£ 7^ 0, M^ ) are shown in 
Figure (jSt). Although the amplitude of the regular wavefunction is negligible compared 
to the s-wave case, the one of the irregular wavefunction is not. This means that the 
contribution of an £ 7^ channel to the shift, which involves the product of these 
wavefunctions (see equation (J23|) ) may not be negligible compared to the s-wave case. 
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The contributions of the S, Ms; i ^ 0, Mi channels are estimated with a very good 
accuracy by neglecting the electronic interaction potential compared to the centrifugal 
potential when integrating the radial Schrodinger equation. This shows that these 
contributions are independent of the interaction potential and thus independent of the 
s-wave scattering length. As the range Rmin < r < Rmax lies in the classically forbidden 
region in the centrifugal barrier, the exact phase shift in the i > wave does not need 
to be known. The regular and irregular wavefunctions are very well approximated by 
regular and irregular spherical Bessel functions jiik'^^r) and ye{k'^igr). 



^5,M,;.^o,M,(r) = J-^ (kTisr) Ji{ktisr) (2^ 



TT^Ms 



^SMs;m,MAr) = \l:^ (kMsr) yiikTisr) (29) 



^Ms 



Using the expansion of the wavefunctions ()28|) ()29|) for small k'^ , we find the 
following dependence on scattering energy and magnetic field : 

/CO pri 



(2£+l) r{ ' 2 (2£-l)(2^+3) 



5.3. Results 



The results (j^ and (IHUj) for each channel S, Ms; i, Mi are summed according to © to 
give the final result. Table Ogives the a dependence of the shift in the case E^o = 0, 
-8 = and a~ polarization of light. It shows the partial contributions of the £ = 0, 2, 4 
collisions, their sum, and the total shift for a = 143 Oq ^21- The s-wave contribution to 
the shift is dominant, and according to section I3?11 it is also independent of the magnetic 
field. As a consequence, these results are quantitatively similar to what is obtained in 
the experimental conditions of the experiment described in pp. 

Table 1. Contributions to the shift in MHz/W.cm"^ due to s,d, and g-wave colhsional 
channels in case of couphng by (t~ hght, with a magnetic field B = 0, and with 
a scattering energy Eoo = 0. a is the scattering length for the ground-state ^S+ 
potential expressed in unit of oq. The last column gives the result for a — 143 ag 
which is the value measured in |12| . 

s-wave d-wave g-wave total total (a — 143 oq) 

v=0 -16.75-1-0.0862 a 

v=l -13.47-1-0.0310 a 

v=2 -38.16-1-0.1001 a 



1.89 


-0.06 


-18.70-H0.0862 a 


-6.37 


2.64 


-0.02 


-16.13-F0.0310a 


-11.70 


5.81 


-0.01 


-43.88-1-0.1001 a 


-29.57 
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Figure 6. Ratio of the shifts 5'"^^/6'"^^ and (J^^^/J"^" as a function of the scattering 
length a in case of a cr^ coupling and with a magnetic field 5 = 0. At a « 210 oq, S""^^ 
changes sign. For a < 210 ag the level v = is shifted towards red detuning which was 
also observed experimentally 0. 



As mentioned earlier, the level shift is a linear function of the coupling elements 
r2| M -i M 1 ^"^^ ^o ^ linear function of the laser intensity /. However, the calibration of the 
laser intensity is difficult to achieve experimentally. The ratio of the shifts for different 
levels is easier to measure accurately. Figure IHl shows the ratio of the shifts 5""^^ /5^^^ and 
^v=2 i^v=Q g^g ^ function of the scattering length a in case of the coupling with a purely 
a~ light and for a magnetic field 5 = 0. The singularity at a ^ 210 ao is associated to 
the cancellation and the change of sign of 5""^^ . The curves of Figure El adapted to the 
B ^ Q case have been used in PP to deduce from these measurements a range of values 
for the scattering length a. The measured ratios in P are 5^^^ /5^^^ = 1.70 ± 0.14 and 
^v=2^^v=o _ 2 2 _|_ g-j?^ g^j^^ correspond to the left part of the figure where all 6"" are 
towards red detuning. Let us notice that Table [T] is given for fc^^^g = 0; so that the 
linear dependence of the s-wave contribution to the shift is linear in a for any value of a. 
In the experiment described in [T], fcMs=2('^~^) << ^ ^^^ Rmin < r < Rmax and a in the 
expected range for spin-polarized metastable helium. Therefore, according to section 
15.11 the s-wave contribution to the shift can still be considered as a linear function of a. 

Table El shows the results obtained with a"*" polarization and i? = as a function 
of the scattering length a. As each ground state channel is coupled either by a" or a"^ 
light, the corresponding contributions to the shift can be weighted by the amount of laser 
intensity of each polarization and then added up to account for a given experimental 
situation where both polarizations are present. We remind that as the experiment 
described in [T] uses spin polarized atoms which form a J = 2,Mj = 2 colliding pair. 
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and since we are interested in a J' = 1, Mji = 1 excited state, the PA laser must 
necessarily contains a a~ component in order to excite the atoms. Non a~ components 
of light do not excite the atoms, but induce an extra shift of the excited states. 

Table 2. Contributions to the shift in MHz/W.cm"^ due to s,d, and g-wave colhsional 
channels in case of couphng by cr+ hght, with a magnetic field B = 0, and with 
a scattering energy Eoc — 0- a is the scattering length for the ground-state ^I]+ 
potential expressed in unit of oq. The uncertainty in the s-wave contribution is 
due to the uncertainty on the scattering length of the ground state ^E+ potential 
fli = 35 ± 15 ao(see section |OJ. The last column gives the results for ai — 35 ao and 
a = 143 ao H!]. 

s-wave d-wave g-wave total total (a = 143 ap) 

v=0 -5.86±0.03-h0.0144 o ^3^9 ^(101 -9.36±0.03+0.0144 a -7.36± 0.03 

v=l -6.89±0.02-h0.0052 a -4.10 -0.004 -10.99±0.02-h0.0052 a -10.25± 0.02 

v=2 -17.80±0.06+0.0167a -8.70 -0.002 -26.50±0.06-|-0.0167 a -24.11± 0.06 

In zero magnetic field, we can notice that the results of the last column in Table |21 
are less than 15 percent different from those of Table ^ However in the experimental 
situation of [T], the magnetic field is non zero, and k^^i^^^Rmax — 1 so that the 
approximation mentioned in section 15.11 breaks down. Numerical calculations of the 
scattering wavefunctions show that the s-wave contribution to the shift is strongly 
suppressed for t> = 2 with cr+ polarization, and the overall shift is significantly reduced 
in absolute value compared to the a~ coupling case. This is illustrated in Figured The 
shifts for f = 0, 1, 2 are plotted versus the polarization of light in the conditions of the 
experiment described in P, where linear polarization was investigated in addition to 
a~ polarization. 

5.4- Additional light shifts 

All the previous calculations deal with the light-induced frequency shift of the vibrational 
level which is excited by the PA laser in the Ojj^ potential. The initial scattering state 
in the ^EJ" potential is also shifted by light. When the two colliding atoms are far 
apart, the non resonant laser excitation, detuned by an amount A < from the atomic 
S — Pq transition, induces a shift of the internal atomic energies which is of the order 
of n^/A, where Vtoo is the Rabi coupling at large distance. For the photoassociated 
bound level w = 0, 1, 2 | A| > 200 MHz is much larger than the width 7/27r = 1.62 MHz 
of the excited atomic P state, so that 7 does not appear in the expression of the shift. 
Since A is negative, the shift is negative and produces an additional blue shift of the 
photoassociation line. This shift is found to be less than 20 percent of the shift calculated 
in section IKTHl It depends on the photoassociated bound-level as it depends on A. We 
calculated the quantity 5"" /{Vt^^/ A'") in the conditions of the experiment described in 
[Ij, where 5^ is the shift of the level v calculated in sectional A^ is the laser detuning 
associated to the excited bound level v. This quantity has a very small dependence on 
the bound level v. As a consequence, the ratio of the shifts corresponding to different 
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Figure 7. Influence of the polarization of light on the shift of v = 0,1,2 for a 
magnetic field corresponding to a Zeeman shift of 20 MHz between the Ms — 2 and 
Ms = channels, a scattering energy Eoo/h = 0.1 MHz (experimental conditions of 
the experiment described in 1 ) and a scattering length a = 143 oq [12j . I„± is the 
laser intensity in the a polarization. For v = 2 the shift is more dependent on the 
polarization than for w = 0, 1. This is due to the fact that the value of the integral 
l|23|l is approximately zero for v — 2 and for the scattering wavefunction for the a'^ 
coupled channel S* = 2, Ms = 0; £ = 0, M^ = at energy (fc^g^o)^/^^- Therefore the 
s-wave contribution to the shift is strongly reduced: |(5^^^| is much smaller for a (t+ 
polarization than for a <7~ one. 



vibrational levels, which is used to determine the scattering length (see Figure IH} is 
calculated to be the same within 3 percent, whether this additional shift is taken 
into account or not. The difference is well below the experimental accuracy of the 
measurements described in [Q and we have neglected it. 

A precise treatment of this effect would require a description of the PA as a 
scattering problem, and to relate the experimental signal to the loss of atoms, i.e to the 
loss of unitarity of the ^-matrix. (1 - ZIs/^m' /',m' P5',m^,£',m;;S',Ms;<?,mJ ) iMj- This loss 
rate exhibits a resonance which is shifted from its position in the zero intensity limit by 
the quantity 5 from equations © and (J23|l . whatever the incoming channel 5, Ms] i, Me 
is |H]. The difficulty which is then encountered is that the scattering channels remain 
coupled at infinite atom separation due to the non resonant excitation of the atomic 
S — Po transition. A change of basis must therefore be introduced in order to use 
states which are uncoupled at infinite separation (dressed states), as they should be 
in a S'-matrix description. Some specific features linked to the proper definition of 
the boundary conditions in photoassociation called dressing effects were investigated 
in [^ Eni EH El ■ However, photoassociative frequency shifts were not studied in this 
context. Treating them would require to take into account more excited channels, for 
example the lu repulsive potential connected to the S — Pq asymptote. It is involved 
in the non-resonant laser excitations at large distance since it is degenerate with the 
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purely long range 0+ potential in the limit of infinite interatomic distance. The present 
accuracy of the experiment does not yet require a more precise calculation of this type. 

6. Conclusion 

In this article we have given a detailed theoretical analysis of the light-induced frequency 
shifts in the PA of "^He* atoms. We have considered the case where a pair of spin 
polarized atoms are excited into a molecular state in the purely long range J = 1,0^ 
potential by absorbing one photon at a frequency red to the 2^5*1 — 2^Pq transition. 
The coupling of this excited state to ground bound levels and scattering states has 
been considered and the physical role played by theses various contributions has been 
discussed. We have also emphasized the importance of the contribution of the £ 7^ 
scattering channels. Explicit expressions of the light shifts have been derived from an 
effective Hamiltonian approach. Results have been given for several vibrational levels 
(f = 0, 1, 2) in the J = 1, 0^ potential and for both a~ polarization of the laser light, 
which is actually absorbed by the spin polarized atoms, and cr"^ polarization, which is 
not absorbed but contributes to the shift through the coupling to various channels. The 
explicit calculation of the shift ratios for two different values of v, which are independent 
of the laser intensity, provided an appropriate tool to interpret previous experimental 
results of our group and allowed us to derive a value for the scattering length a[T]. 

Actually the ^He* atom can be regarded as a model system for simply solving the 
problem of the light shifts in a photoassociation experiment. First, the entrance channel 
for the pair of incoming atoms is well defined, as the atoms are spin polarized and the 
temperature is ultra low. Second, the helium atom has no hyperfine structure, which 
limits the number of coupling channels to be considered. In addition, the fact that the 
excited molecular state 0^ is purely long range allows an exact calculation of the excited 
molecular potential and wave functions, and the use of an asymptotic expansion of the 
wave function for the incoming pair of ground state atoms thus avoiding its complexity 
at short interatomic distance. This is why one simply finds a linear dependence of 
the shifts as a function of the s-wave scattering length a. Measuring the light-induced 
frequency shifts is thus a very appropriate method to deduce a value for a without having 
to precisely know the ground state potential. In that respect the light-shift method is 
more direct yet much less precise than the two-photon Raman spectroscopy method 
|12j . in which the relation between the measured energy of the least bound state and 
the value of a has to be infered from the ground state potential. 

Alkaline-earth atoms are other model systems in which light shifts could be studied. 
They are currently the focus of intense studies. They appear as simpler systems than 
alkalis since they have isotopes lacking hyperfine structure. From light shifts one could 
infer complementary or more accurate results on scattering lengths than those given by 
other methods. One could also test photoassociation theories on these simple systems. 
This could for instance apply to Ca, Sr or Yb atoms. The isotopes ^'^Ca, ^^Sr, ^^Sr 
and ^"^^Yb have been studied in j38 | IHn t 1 ^ W H, I42j . The case of ^^Sr is of special interest 
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as it has long range potentials which were investigated in |13|; and the case of ^'^^Yh is 
specific as it displays atomic levels without fine nor hyperfine structures. 
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Appendix A. Symmetrized states for collisions of metastable atoms 

The physical states to be considered for describing the collisions between '^He* atoms 
have to be symmetric under the permutation of the bosonic identical atoms. These are 
obtained using the symmetrization operator (1 + Pah)l1 acting on the non symmetrized 
basis (JD). Using the transformation in the spin space from the |(S'a, S'fe)S', M5) to the 
\Sa-, Ms^] Sh, Ms^) basis, where Sa^S^, and Ms^^Ms^ are the spins of atoms a, h and their 
projection on the lab-fixed axis respectively. Using also a Clebsch-Gordan identity, we 
get : 

\{Sa.S,)S,Ms) = Yl {Sa,Ms^;S,,MsJ{Sa,S,)S,Ms)\Sa,Ms^;S,,Ms,) 

= Yl {-l)'^^''-'{St,Ms,;Sa,MsJ{S,,Sa)S,Ms) 

\Sa,Ms^;Sb,Ms,) 
We deduce the action of Pat on the spin part: 

Pab\iSa,S,)S,Ms) = {-lf'^+''^-'\iS,,Sa)S,Ms) 
The action of Pat on the spatial part is given by: 

Pak\i,Me) = {-lY\i,M,) 

We get the properly symmetrized basis : 

1 + Pk 1 (_l)Sa+Sb-S+i 

-I-^\(Sa,S,)S,Ms;i,M,) = -\{Sa,S,)S,Ms;l,Mi)+^-^ \{S,,Sa)S,Ms;i,M,) 

(A.l) 

In the case of ^ife*, Sa = Sb = I, which implies S' = 0, 1, 2 (singlet, triplet and quintet 
states). We will therefore use a simplified notation: 

\S,Ms-,e,Me) = ^\{Sa,Sk)S,Ms-,l,Mi) + ^-^^^^^-^\{Sb,Sa)S,Ms-,e,M,) (A.2) 

The formula ()A.2|) also shows that singlet (5* = 0) and quintet (5* = 2) states collide in 
even i partial waves, and triplet states in odd i partial waves. 

We now investigate the properties of the collisional state under the inversion of the 
electrons coordinates through the center of charge of the molecule (gerade/ungerade 
symmetry). In the case of homonuclear dimers, the center of charge and the center of 
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mass are the same. As both the ^He* atom and its nucleus are bosons, the exchange of 
the atoms and the exchange of the nuclei of the quasi-molecule should not change the 
wavefunction of the pair of atoms. The product of both operations is equivalent to the 
exchange of the electronic clouds and should therefore not change the total wavefunction. 

In the space coordinates the exchange of the electronic clouds is the product of the 
inversion of the electronic coordinates through the center of mass of the corresponding 
atom, and the inversion of the electronic coordinates through the center of mass of the 
molecule (gerade/ungerade symmetry). As the atoms are in a 2^ Si state with an even 
spatial electronic wavefunction, the first operation is the identity. The second operation 
multiplies the wavefunction by (—1)'", where w = for a gerade state, and w = 1 for 
an ungerade state. 

In the spin space, the exchange of the electronic clouds multiplies the wavefunction 
by (-l)^''+^'>-s = (-1)^. 

As a consequence, (—1)"'+'^ = +1 due to the bosonic nature of the '^He* atom and 
of its nucleus. Singlet (5* = 0) and quintet {S = 2) states are gerade states {w = 0) 
whereas triplet states {S = 1) are ungerade states {w = 1) 

Appendix B. Light-induced couplings 

In this appendix we calculate the radiative coupling operators Qegi'f') introduced in 
section IH and connecting the excited channel e : (0^; J' = 1, Mji = 1) to the ground 
state channels g : (5, M^; £, Mi) given in equation (0) for a a^ excitation and in 
equation (jH)) for a o""*" excitation. More precisely, the ^egij) are obtained by taking the 
matrix elements of the laser interaction Hamiltonian in the electric dipole approximation 



27r//c d.e between (e| and \g). 

^eg{r) = ^sms-AmA^) = -v/2^(J' = l,Mj, = l;0+|d.e|5, M5;£,M,) (B.l) 

In flB.l|) . I is the laser intensity, c is the speed of light, d the dipole operator of the 
quasi-molecule, e is the polarization vector of light. Note that the matrix element ()B.1|1 
is taken over the electronic variables and the angular variables of the molecular axis, so 
that ^sMs'/mX''^) remains an operator for the one-dimensional variable r, the distance 
between the two nuclei. Using the decomposition (jHl) of the excited state, ^s,Ms-/Mei''^') 
can be expressed as : 

4 

^sMs;ImX^) = -\/27r//c^a,(r)(J' = l,Mj, = l,^|d.e|5, Ms; /, M;) (B.2) 

In equation ()B.2j) . the index i stands for one of the electronic states ^S^,^nu,'^n„,^S+ 
appearing in the expansion of the excited molecular state \J' = l,Mj/ = 1,0^), the a, 
coefficients are found by the diagonalization of the electronic Hamiltonian for a S + P 
molecule (see subsection 13. 3^ . The (J' = l,Mj/ = l,i\d.e\S, Ms; I, Mi) matrix elements 
are needed to calculate (IB. 21). 
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The electric dipole operator acts only on the orbital angular momentum of the 
electrons, and can therefore be conveniently expressed in the Hund's case (a) basis for 
the rotating molecule ()B.3jl . in which its projection on the molecular axis is well defined. 



2 7 + 1 

\J,Mj,n,S,J:,A,w) ^ ^ ^^Di,*^^{a,(3,^ = 0)\S,J:,A,w) (B.3) 

J = L + S + i is the total angular momentum, which is the sum of the orbital angular 
momentum of the electrons L, the spin of the electrons S and the angular momentum 
for the rotation of the atoms i which has a zero projection on the molecular axis. Mj 
is the projection of J on the lab-fixed quantization axis, A and E the projection of L 
and 5* on the molecular axis, fi = A + S is the projection of J on the molecular axis. 
The transformation of the wavefunction under the inversion of the electronic space 
coordinates is characterized by a quantum number w such that (—1)"' = 1 (gerade 
symmetry) or (—1)"' = —1 (ungerade symmetry). Df^* ^ is a matrix element of the 
(2 J + 1) X (2 J + 1) unitary rotation matrix which is a function of the Euler angles 
Q!,/3,7 and which defines the relationship of the space-fixed lab frame with that of the 
body-fixed molecule frame. 7 = is a convention for linear molecules for which two 
angles are enough to describe the rotation. 

The states wich are involved in the decomposition ()B.2|) are related to the basis 
(|R3|l by the relations (JR4|l : 

|J' = l,Mj, = 1,5 S+): 11,1,0,2,0, 0,1) 

11,1,0,0,0,0,1) 

1/v^ (|1, 1, 0, 1, 1, -1, 1) - |1, 1, 0, 1, -1, 1, 1)) 

I/V2 (|1, 1, 0, 2, 1, -1, 1) + |1, 1, 0, 2, -1, 1, 1)0.4) 

In flB.4|) . the linear superposition for the n„ states is formed so that the electronic states 
have (+) symmetry under the reflexion of the space and spin electronic coordinates about 
a plane containing the molecular axis, and are thus part of the 0^ subspace. 

The dipole operator d is best expressed in the body-fixed molecular frame, and e 
in the space-fixed frame. The product d.e can be expressed using a rotation matrix D, 
using the quantum numbers n and m which refer to the projection of d on the molecular 
basis and e on the space fixed basis (m = 0, +1, —1 for 7r,o"+,o"~ polarization) |44j . 

d.e= Yl i-ire.^dr,Dl\ia,P,0) (B.5) 

m,n=0,itl 

The projections in the molecular frame dn=o,±i of the operator d act on the electronic 
part |S', S,A = 0,w) of a Hund's case (a) ground state according to equation ()B.6|) . 
where dat is the atomic electric dipole moment for the S — P transition [44J. 

{S',E',A',w'\dn\S,E,A = 0,w)= '-±t:l^r:^dat6s,s'Sj:,j,'6A>,n (B.6) 

Using ()B.3|) . ()B.5|) . ()B.6|1 . we can form the expression of the action of the electric dipole 
operator in the Hund's case (a) basis for the rotating molecule. 



|j' = l,Mj, = 1,3h,) 
|j' = l,Mj, = 1,5H„) 



{J',M'j,n',A',S',J:',w'\d.e\J,Mj,Q,A = 0,S,J:,w)- ./(2-^+l)(2-^' + l) 



167r2 
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l + i-l 



iw' — W + l 



2-K 



—E. dat5s,S'^T.,T.'^K\n 2_^ ^-m{ — ^ 

* m,n 

/•TT 

da / sinl3dl3Di,^ ^,ia, /3,0)Dl\ia, (3,0)Di,; f,(a, ^0) 



Uo Jo 

The evaluation of the integral over rotation matrix elements gives |44j : 



(J', M'j, n', A', S', S', w'\d.e\J, Mj, f], A = 0, S, S, w) 
J2^-mhi)"'{J,Mj-l,m\J',Mj,){J,Q;l,n\J',Q') 

m,n 



2J+1 
2J' + 1 



v^ 



-'^at'^5,S"^E,S"5A',n 



(B.7) 



Now, the transformation of the electronic ground state between the basis 
\S,Ms]£,Me) (fXT2|l and the Hund's case (a) basis (JR8|l with A = is needed. The 
rotation in the spin states between the lab-fixed frame and the molecular frame is 
expressed in the relation ()B.8|) . 

1^; Ms) = J2 D'm^Ac^, a 0) |5, S) (B.8) 

A state \i, Me) is associated to the spherical harmomic Y^^Mi which can be expressed 
using a rotation matrix element (Appendix C in fl^)'. 



|£;M,)^F,,M,(A«) 



2^ + 1 p£. 

Air 



D%^,{a,l3,0) 



(B.9) 



In the electronic ground state, the orbital angular momentum of the electrons is zero 
so that A = 0. The property under gerade symmetry is given by (—1)"'+'^ = +1 (see 
appendix A). We write therefore using ()B.8j) and (JR 



\S,Ms;i,M,) 



2i + l 
A-n 



^koKAO) 5^Z}f,;,s(«,A0)|S,S,A = 0,^ = S) (B.IO) 



The following scalar product is deduced 

(J,Mj,f],^',S,A = 0,M;|^,M5;£,Mf) = 



(2J+l)(2f+l)(l + (-l)'"'+^) 
167r2 2 



5,S' 



27r 



(B.ll) 



da / ^in(3d(3Di^^^{a,(3,Q)Dil,{a,(3,Q)Dl,l^{a,(3,Q) 

'0 ^0 

Performing the integration of the product of rotation matrix elements, we obtain the 
transformation ()B.12|) between the basis ()A.2|) and ()B.3|) : 



(J, Mj, n, S', ^,A = 0,w\S, Ms; i, Me) 



2i + l 
2J+1 



Ss,s'S^,s{S,Ms;i,Me\J,Mj){S,^;i,0\J,Q) 

(B.12) 
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Using ()B.12|) and ()B.7|) . it is now possible to calculate the action of the electric 
dipole operator on the \S, Ms; i; Mi) projected in the Hund's case (a) basis ()B.3|) 

{J',M'j,Q',A',S',J:',w'\d.e\S,Ms;i,M,) = Y, Yl {-^T^-m^d,t ^ ^ 

m.n J,Mj,n,T, 



6A>,Js,s'Sj,,^'^^jr^{J,Mj-,l,m\J',M'j){J,n-,l,n\J',n'){S,Ms-,eMe\J,Mj) 

The summation can be simplified using the following remarks: 

• The fourth Clebsch-Gordan coefficient is non zero for fi = S. Because of the 
coefficient ^s^s', we have fl = T,' 

• The third Clebsch-Gordan coefficient is non zero for Mj = Ms + Mi 

• The second Clebsch-Gordan coefficient is non zero for n = Q' — Q = Q' — Ti' = A' 

• The first Clebsch-Gordan coefficient is non zero for m = Mji — Mj = Mji — Ms — Me 

This gives finally equation ()B.13|) 

( /, M'j, n', A', S', S', w'\d.e\S, Ms; i, Me) = Y,i-^)''''''''~^''^-M,,+Ms+M,V2da 



''at 



I A_ (r]w'+s-i\ loP+l 

{J,Ms + Me;l,Mr-Ms-Me\J\Mj.) 



2 y V 2 J' + 1 

(J, S'; 1, AV', ^'){S, Ms; IMfp, Ms + Me){S, S'; I, 0| J, S') (B.13) 

As a conclusion, using ()B.13|) . and the relations ()B.4|) . it is possible to calculate 
(JR2|l . The I J' = 1, Mj> = 1,^ n„) state does not contribute to the sum (JR2|l due to the 
selection rule on the u/g symmetry ()B.6|) and the fact that ground electronic states are 
gerade states (see [Appendix A\) . 
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